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ABSTRACT 

Several transiting super-Earths are expected to be discovered in the coming few years. While 
tools to model the interior structure of transiting planets exist, inferences about the composition 
are fraught with ambiguities. We present a framework to quantify how much we can robustly infer 
about super-Earth and Neptune-size exoplanet interiors from radius and mass measurements. We 
introduce quaternary diagrams to illustrate the range of possible interior compositions for planets 
with four layers (iron core, silicate mantles, water layers, and H/He envelopes). We apply our model 
to CoRoT-7b, GJ 436b, and HAT-P-llb. Interpretation of planets with H/He envelopes is limited 
by the model uncertainty in the interior temperature, while for CoRoT-7b observational uncertainties 
dominate. We further find that our planet interior model sharpens the observational constraints on 
CoRoT-7b's mass and radius, assuming the planet does not contain significant amounts of water or gas. 
We show that the strength of the limits that can be placed on a super-Earth's composition depends 
on the planet's density; for similar observational uncertainties, high-density super-Mercuries allow 
the tightest composition constraints. Finally, we describe how techniques from Bayesian statistics 
can be used to take into account in a formal way the combined contributions of both theoretical 
and observational uncertainties to ambiguities in a planet's interior composition. On the whole, 
with only a mass and radius measurement an exact interior composition cannot be inferred for an 
exoplanet because the problem is highly underconstrained. Detailed quantitative ranges of plausible 
compositions, however, can be found. 

Subject headings: planets and satellites: general, planetary systems, stars: individual (CoRoT-7, 
GJ 581, GJ 436, HAT-P-11) 



1. INTRODUCTION 

Over two dozen low-mass exoplanets with masses less 
than 30 Earth masses are knowrfj. As their numbers in- 
crease, so does the probability to uncover a population 
of transiting low-mass exoplanets. The fir st transiting 
super -Earth exoplanet has been discovered (|Leger et al.l 
2009)-based on the young history of exoplanets once 
one example of new type of object is discovered many 
more soon follow. Now that we are on the verge of dis- 
covering a good number of low-mass transiting planets 
(iBaglin et ai ; l200ftjBorucki et al.l2008l;llrwin et ; al.l2009t 
iMavor et al.ll2009bt lLovis et al.ll2009D . methods to con- 
strain their interior composition from observations are 
required. 

A good example of why quantitative methods to 
constrain pl anetary interior compositions are nee ded 
is G J 436b (IButler et al.1 l200i IGillon et al.l l2007bT). a 
Neptu ne-mass (M p = 23.17 ± 0.79M e ,; iTorres et al.l 
2111), Neptune-size (R p = 4.22t£;^i? e ; ITorres et all 



(pOOl ) planet in a 2.6-day period around an M2.5 star. 
Initially, because of its similarity to the physical pro- 
portions of Neptune, IGillon et al l (|2007bl ) assumed that 
GJ 436b was composed mostly of ices. Others showed 
that it could instead be compose d of a rocky interior with 
a more massive H/He envelope (jAdams et al.ll2008h . 
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Previously, iValencia et al.l ()2007al ) introduced ternary 
diagrams to constrain the int erior composition of su per- 
Earths without gas envelopes. iZeng fc Seagerl ()2008l ) pre- 
sented a detailed description of the functional form of 
the ternary diagram interior composition curves. Super- 
Earths are loosely defined as planets with masses be- 
tween 1 and 10 Earth masses that are composed of 
rocky or iron material. While the terms "mini-Neptune" 
or "Neptune-like" are not in common usage, they re- 
fer to planets with significant gas envelopes. Others 
have modeled evolution of Neptune-mass planets to pre - 
dict radii (e.g. iFortnev et all 120071 : iBaraffe et all 120061 ). 
iFigueira et al.l ()2009f ) used planet formation and migra- 
tion models to suggest interior compositions for GJ 436b. 

In this paper, we aim to quantify the constraints placed 
on a low-mass exoplanet 's interior structure by transit 
and radial velocity observations. We use a planetary 
structure model to explore the range of plausible inte- 
rior compositions that are consistent with a given pair 
of mass and radius measurements, independent of planet 
formation arguments. We extend previous work by in- 
cluding the possibility of a gas envelope and by consid- 
ering a range of mantle iron enrichments. Our model 
of low-mass planet interiors includes an iron core, sili- 
cate mantle, water ice layer, and H/He layer. To plot 
four-layer interior compositions we introduce quaternary 
diagrams, an expansion of ternary diagrams into three 
dimensions. Finally, we present a new framework to com- 
bine both model and observational uncertainties in a rig- 
orous way using Bayesian techniques when interpreting 
the interior composition of a transiting exoplanet. Our 
overall goal is to be able to interpret planetary mass and 
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radius observations with a quantitative understanding of 
the effects of model uncertainties, observational uncer- 
tainties, and the inherent degeneracy originating from 
the fact that planets of differing compositions can have 
identical masses and radii. 

We describe our planetary interior structure model in 
§ We introduce quaternary diagrams in § [31 In § [H we 
apply our model to constrain the compositions of low- 
mass exoplanets. In § [5j we describe how Bayesian tech- 
niques may be applied to the problem of drawing infer- 
ences about an exoplanet's interior from measurements 
of the planet's mass and radius. Discussion and conclu- 
sions follow in § |5] and [7J 

2. MODEL 

2.1. Model Overview 

We consider a spherically symmetric differentiated 
planet in hydrostatic equilibrium. With these assump- 
tions, the radius r (to) and pressure P(m), viewed as 
functions of the interior mass to, obey the coupled dif- 
ferential equations 
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where p is the density and G is the gravitational con- 
stant. Equation (JT]) is derived from the mass of a spher- 
ical shell, while Equation @ describes the condition for 
hydrostatic equilibrium. The equation of state (EOS) of 
the material 



P = f(P,T) 



(3) 



relates the density p (m) to the pressure P (m) and tem- 
perature T (to) within a layer. We allow our model plan- 
ets to have several distinct chemical layers ordered such 
that the density p (to) is monotonically decreasing as to 
increases toward the planet surface. Throughout the rest 
of this work, we shall use Xi — Mi /M p to denote the frac- 
tion of a planet's total mass M p in the ith layer from the 
planet center (i = 1 denotes the innermost layer). 

To model a planet having mass M p , radius R p , and 
a specified composition {a^}, we employ a fourth-order 
Runge-Kutta routine to numerically integrate Equa- 
tions ((TJ) and © for r (to) and P (to) from the outer 
boundary of the planet (to = M p ) toward the planet cen- 
ter (to = 0). We describe our sche me for setting the exte- 
rior boundary conditions in § 12.31 We impose that both 
P and r are continuous across layer boundaries. At each 
step in the integration, the EOSs and temperature pro- 
files described in § 12.21 are used to evaluate p (to) . 

The planet parameters {M p , R pi {xi}} in fact form an 
overdetermined system; there is a single radius R p that 
is consistent with M p and {xi\. For a given mass and 
composition, we use a bisection root-finding algorithm 
to iteratively solve for the planet radius Rp that yields 
r (m = 0) = upon integrating Equations §B) and ([!} to 
the planet center. We stop the iteration once we have 
found R p to within 100 m. Alternatively, in some appli- 
cations it is convenient to be able to stipulate the planet 
radius (for instance when exploring the range of compo- 



sitions {x^ allowed for a confirmed transiting exoplanet 
of measured mass and radius). In these cases, we use 
a bisection root-finding algorithm to iteratively solve for 
the mass ratio of the inner two material layers (xz/xi) of 
the planet given M p , R p , and valid mass distribution in 
the outer layers of the planet {xi | i > 2}. We stop this 
iteration once we have found x\ and X2 to within 10 -10 . 

We increase the achievable accuracy in the composi- 
tion of our modeled planets and the stability of this it- 
erative process by employing the Lagrangian form of the 
equations of structure. With mass to as the indepen- 
dent integration parameter, we can take a partial mass 
step at the conclusion of each layer i to ensure that the 
specified value of Xi is precisely obtained. Within each 
layer, we employ an adaptive mass step-size such that 
each integration step corresponds to a radius increment 
of approximately 100 m. An adaptive step-size is neces- 
sary because both Equations ([IJ and ([2]) diverge as r — >• 
and to — y 0. 

2.2. Material EOS and Thermal Profiles 

In this section, we describe the EOS and thermal pro- 
file T (to) assumed for each material layer. 

We allow for the presence of an outer gas envelope in 
our modeled planets. We use t he H/He EOS with h elium 
mass fractio n Y = 0.28 from[Saumon et al.l (|1995f ). As 
mentioned in Ad ams et al.l we ignore the "plasma 

phase transition" in the H/He EOS. To set the thermal 
profile we divide the H/He layer into two regimes: a thin 
outer radiative layer and an inner convective layer. 

In the radiative regim e of the gas e nvelope, we em- 
ploy the analytic work of Hansen ^2 008Tl t o approximate 
the temperature profile. iHansenl (|2008t ) considered a 
plane-parallel atmosphere in radiative equilibrium that 
is releasing heat flux generated in the planet interior 
while also being irradiated by a mono-directional beam 
of starlight. He solved the gray equations of radiative 
transfer with a 'two-stream' approach, allowing the in- 
coming optical stellar photons to have a different opacity 
and optical depth than the infrared photons reradiated 
by the planet, and obtained a temperature profile 
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In the above equation, T is the atmospheric temperature, 
r is the infrared optical depth, 7 is the ratio between the 
optical and infrared optical depths, po is angle cosine 
of the incoming beam of starlight relative to the local 
surface normal, T c g is the effective temperature of the 
planet in the absence of stellar irradiation, and To char- 
acterizes the magnitude of the stellar flux at the orbital 

distance of the planet ^F* (R^/a) 2 = crTo 4 ^ . While po 

varies over the surface of the planet, our planet model is 
one-dimensional spherically symmetric model. We adopt 
a single fiducial value of ^0 = 1/2 (the average of po over 
the day hemisphere) when calculating the temperature 
profile of the radiative gas layer. Equation ((4]) yields the 
temperature in the radiative gas layer as a function of the 
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(infrared) optical depth. The variation of optical depth 
r, with interior mass m obeys 
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where k is the opacity. In the radiative regime of the 
gas layer, we integrate Equation ([5]) along with Equa- 
tions (QJ and ([2]). For k, we use tabulated Rosseland 
mean opacities of H/H e at solar abundanc e metallicity 
([M/H] = 0.0) from lFreedman et all (f2008h . 

In our model gas layer, we allow for the presence of an 
inner adiabatic regime within which energy transport is 
dominated by convection. Neglecting the effects of con- 
duction and diffusion, we take the temperature profile 
in the convective layer to follow the adiabat fixed to the 
specific entropy at the base of the radiative regime. The 
transition between the radiative and convective regimes 
is determined by the onset of convective instabilities. An 
adiabatically displaced fluid element in the gas layer will 
experience a buoyancy force tending to increase its dis- 
placement if 
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where the density p = p (P, s) is viewed as a function 
of the pressure P and specific entropy per unit mass s. 
Whenever Equation © is satisfied, the gas layer is unsta- 
ble to convection. In the H/He EOS from lSaumon et aD 
(1995), the adiabatic gradient (dT/dP) s is positive for 
all values of P, T, and He mass fraction Y. It thus suf- 
fices to test for ds/dm < to define the outer boundary 
of the convection regime. As we integrate Equations ([T]), 
([2]l. and ([5]) from the planet exterior inward, we transi- 
tion from the radiative regime to the convective regime 
once ds/dm < 0. 

In the interior solid layers of the planet, we neglect the 
temperature dependence of the EOS. Thermal effects in 
the solid laye rs of a planet hav e a small effect on the 
planet radius (jSeager et al.ll2007l ) justifying the assump- 
tion of a simplified isothermal temperature profile. For 
every solid material considered in this study (Fe, FeS, 
Mgi-yFeySiOa, and H2O) we use the EOS data sets from 
ISeager et al.l ([2007D derived by combining experimental 
data at P < 200 GPa with the theoretical Thomas- 
Fermi-Dirac EOS at high pressures, P > 10 4 GPa. 

2.3. Exterior Boundary Condition 

In our model (described in § I2.1[) . the exterior bound- 
ary of the planet sets the initial conditions for integrat- 
ing the equations of structure. In the absence of a gas 
layer, we take the pressure to be at the solid surface of 
the planet (m = M p , r = R p , P = 0). For planets hav- 
ing gas layers, we use a simplified constant scale height 
atmospheric model to choose appropriate exterior bound- 
ary conditions on the pressure P and optical depth r at 
r (M p ) — R p as elaborated below. 

To physically motivate our choice of exterior boundary 
conditions for planets with gas layers, we make several 
simplifying approximations about the properties of the 
planet gas layer in the neighborhood of the measured 
planet radius R p . We assume that in this region the gas 
layer can be approximated as an ideal gas, so that 



P = 



peS 



(7) 



where p c s is the effective molecular mass of the gas. We 
further assume that the outer atmosphere of the planet is 
isothermal. This i s consist e nt wit h the radiative temper- 
ature profile from [Hansen (2008) (Equation (j4])), which 
is largely isothermal for r<l. We also neglect varia- 
tions in the surface gravity g = GM / R 2 over the range of 
radii being considered. Finally, to account for the pres- 
sure dependence of the opacity, we assume a power-law 
dependence 
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where logC = -7.32, a = 0.68, and = 0.45 a re de- 
termined by fitting to the iFreedman et al.l (|2008|) tabu- 
lated opacities (with all quantities in SI units). These 
assumptions, when coupled with the equation of hydro- 
static equilibrium (dP/dr = —pg) and the definition of 
the radial optical depth (dr/dr = —up), yield an expo- 
nential dependence of both P and r on r, 



p( r ) = p Re -(r--Rp)/Kp j 
r ( r ) =rBe -(a+l)(r-flp)/*P j 

with the pressure scale height Hp given by 
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and the pressure and optical depth at R p (Pr and tr, 
respectively) related by 
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It is important to maintain a direct correspondence 
to observations when defining the radius of a gas-laden 
planet in our model. Planet radii are measured observa- 
tionally from transit depths and thus reflect the effective 
occulting area of the planet disk. We denote the optical 
depth for absorption of starlight through the limb of the 
planet r t (y), where y is the cylindrical radius from the 
line of sight to the planet center. In our models we define 
the transit radius R„ to occur at 



Tt(Rp) = 1. 



(13) 



We use a development similar to that in Hansen (2008) 
to relate the transverse optical depth through the limb 
to the radial optical depth r. Integrating along the line 
of sight through the planet limb, the transverse optical 
depth for starlight is given by: 



, n(r ) p (r) 
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The right-hand side of Equation (fl4"|) is obtained by rec- 
ognizing that for y ~ R p ;§> Hp/ (a + 1) only values of r 
with (r — y) ^ y contribute significantly to the integral 
due to the exponential decay of the integrand. We ob- 
tain exterior boundary condition on r by combining our 
model definition of the transit radius (Equation. (|13[) ) 
with Equation (Ti"4"|) : 
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The boundary condition for pressure follows from tr us- 
ing Equation (1T21) . 



2.4. Model Parameter Space 

In this section, we describe our procedure for choosing 
value ranges for 7, To, and T e s that describe the atmo- 
spheric absorption, stellar insolation, and the intrinsic 
luminosity of exoplanets simulated with our model. 

The parameter 7 in Equation Q denotes the ratio of 
the gas layer's optical depth to incident starlight over its 
optical depth to thermal radiation. At large values of 7 
the starlight is absorbed high in the atmosphere, while 
at small values of 7 the stellar energy penetrates deeper 
into the atmosphere. We adopt a fiducial value of 7 = 1, 
but also consider values spanning from 0.1 to 10. In this 
way, we encompass a wide range of possible absorptive 
properties in our model H/He envelopes. 

In Equation (|4]), /zooTo represents the stellar energy 
flux absorbed (and reradiated) locally at a given point on 
the planet's irradiated hemisphere. The stellar insolation 
impinging on a planet can be calculated with knowledge 
of the host star's luminosity L* or spectral class, and of 
the scmimajor axis a of the planet's orbit. The fraction 
of this energy that is reflected by the planet and how the 
energy that does get absorbed is distributed around the 
planet's surface area, however, remain unknown for the 
super-Earth and hot Neptune planets considered in this 
paper. Our parameterization of the energy received at 
the planet from the star is further complicated by the 
fact that we are using a spherically symmetric planetary 
model, whereas the effect of stellar insolation varies over 
the planet surface. We take these uncertainties into ac- 
count by considering a range of plausible Tq values for 
each planet. For our fiducial value, we use the equilib- 
rium temperature of the planet assuming full redistribu- 
tion and neglecting reflection 
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Similar fiducial choices of To have been made in 
other studies that used Equation ^ to describe the 
gas layer temperat u re profiles of low -mass exoplanets 
ij Adams et all l200l iMiller-Ricci et all 12001 ) . By con- 
sidering reflection of starlight by the planet in addition 
to full redistribution, we set a lower bound on Tq: 



Tq = 



L, (1 - A) 



1/4 



(17) 



where A is the planet's Bond albedo. A plausible upper 
limit of A = 0.35 is chosen; all the solar system planets 



except Venus have Bond albedos below this value. Fi- 
nally, to establish an upper limit on To we neglect both 
redistribution and reflection and take 
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(18) 



This upper b ound co r respo nds to the formal definition 
of T used by lHansenl (|2008D to derive Equation @. 

A planet's intrinsic luminosity (produced by radiogenic 
heating and by contraction and cooling after formation) 
is another important component of the planetary energy 
budgets. In Equation (U), T c ft parameterizes the heat 
flux from the planet interior entering the planet gas layer 
from below, Fi nt — crT c s A . Within the plane-parallel 
gas layer assumption, we can relate T e g to the intrinsic 
luminosity, Li nt , of the planet 
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We require a scheme to constrain the intrinsic luminosi- 
ties of low-mass exoplanets. 

A full evolution calculation, modeling the energy out- 
put of a planet as it ages after formation, is outside of the 
scope of this work. There are many physical effects (in- 
cluding phase separation, chemical differentiation, chem- 
ical inhomogeneities, irradiation, radiogenic heating, im- 
pacts, geological activity, tidal heating, and evaporation) 
that can influence the thermal evolution of a planet and 
flummox attempts to predict a planet's intrinsic lumi- 
nosity (see § 16.31 for a full discussion) . Additionally, the 
ages of the planet-hosting stars considered here (and of 
the planets that surround them) are very poorly con- 
strained. This severely limits the insights that a cooling 
simulation could yield into the planets' intrinsic lumi- 
nosities. Instead of directly simulating planetary evolu- 
tion, we take an approximate scaling approach to bracket 
plausible values for the intrinsic luminosities of low-mass 
exoplanets. 

We use pl anet evolution tracks modeled by 
IBaraffe et all (12001 to constrain the intrinsic lu- 
minosit ies of the gas - laden planets considered in this 
work. IBaraffe et al.l ((2008) modeled the evolution of 
planets ranging from 10 M® to 10 M^, having heavy 
metal enrichments of Z = 2%, 10%, 50%, and 90%, and 
that were either receiving negligible stellar irradiation 
or suffering insolation equivalent to that from a sun at 
0.045 AU. We limit our consideration to the simulated 
irradiated planets that are at least 1 Gyr old and that 
are no more than 1 Mo_. We then fit the intrinsic 
luminosities of this sub-sample of IBaraffe et al.l (2008) 
models to a simple power law in planetary mass, radius, 
and age (t p ): 



log ( — I =ai +8u log [ — — I + bp log ( rjr~ 
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(20) 



The values obtained for the coefficients and their 95% 
confidence intervals are a\ — —12.46 ± 0.05, au„ — 
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1.74 ±0.03, a Rp = -0.94 ±0.09, and a tp = -1.04 ±0.04. 
The fit had R 2 = 0.978 and rms residuals of 0.14 in 
log {Lint/ Lq). For a given planet, we use the measured 
planetary mass, planetary radius, and host star age (a 
proxy for the planet age) with the best-fit coefficients in 
Equation (|20| to calculate a fiducial value for L int . We 
then employ the uncertainties in the fit coefficients, the 
rms residuals, and the range of possible planet ages to 
establish a nominal range of intrinsic luminosities Li nt to 
consider when constraining the interior compositions of 
planets with gas layers. The poorly constrained planet 
age dominates the other sources of uncertainties in its 
contribution to the range of Li nt for all the planets we 
consider. 

Additional limitations on T e ff can be required if the 
nominal range of Li nt determined by the procedure 
above is too broad. At very low values of T e ff (low in- 
trinsic luminosities) the gas layer P-T profile can en- 
ter an unphysical high-pressure low-temperature regime 
(P > 2.5 x 10 10 Pa, T < 3500 K). These conditions, 
under which hydrogen may form a Coulomb lattice or 
a molecular solid, are no t included in the coverage of 
the ISaumon et all (fl995h hydrogen EOS. If necessary, 
we truncate the lower range of T c g values that we con- 
sider to avoid exceed ing the range of applicability of the 
ISaumon et al.l (|1995t l EOS. Out of all the planets consid- 
ered in the work, such a reduction in the range of T c g 
was only required for HAT-P-llb (§l4.4p. 

Adopting a simple scaling approach to estimate T e fj 
allows us to consider a wider variety of possible inte- 
rior compositions than we could by simulating full evo- 
lution tracks. Nonetheless, our procedure to constrain 
T e ff is very approximate. It estimates the intrinsic lumi- 
nosity of a planet from its mass, radius, and age alone. 
The effects of interior composition and stellar irradia- 
tion on a planet's evolution are not addressed. For in- 
stance, because solar system planets are less strongly ir- 
radiated than the transiting planets considered in this 
work, the scaling relation systematically overestimates 
their i ntrinsic luminosities . Further, the extrapolation 
of the iBaraffe et al.l (|2008| ) models to super-Earth-sized 
planets is very uncertain. Although phenomenological, 
the procedure described above provides a consistent way 
to estimate a plausible range of intrinsic luminosities in 
which the span of the range reflects the uncertainties in 
the planet age and thermal history. 

2.5. Model validation 

We have validated our planet interior model by com- 
paring our results with Earth and other models. Our 
fiducial Earth-planet composition is one with a 32.6% 
by mass core consisting of FeS (90% iron and 10% sul- 
fur by mass) and a 67.4% by mass mantle consisting of 
Mg .gFeo.iSi03. For this composition, our model gives a 
radius of 6241 km for a 1 planet. This radius value 
is within 2.2% of Earth's true radius, well within ex- 
pected observational uncertainties for future discovered 
Earth-mass, Earth-sized exoplanets. More importantly, 
our solid planet mod els are not intende d to be accurate 
for such low masses ()Seager et al.ll2007h . because we ig- 
nore thermal pressure. This approximation is much more 
appropriate for more massive planets, where a larger frac- 
tion of the planet's material is at high pressure where 



thermal effects are small. 

We further com pared our mod e l outpu t with the re- 
sults presented in IValencia et all (|2007bf ) . Specifically, 
we reproduced the values in their Table 3 for GJ 876d's 
radius under various assumed bulk compositions . We 
found that for solid planets co mposed of iron and silica tes 
our radii matched those from IValencia et ail (|2007bt ) to 
0.2%. For planets that also included a water layer, our 
radii were within 1%. This is a very reasonable agree- 
ment. The larger discrepancy in the water planet radii as 
compared to the dry-pla net radii stems from differences 
in the EOS for water. See lSeager et al.l (|2007l ) for our cal- 
culations on the water EOS, and a detailed description 
of our EOS choices. 

We tested our model of planets with signific ant gas en- 
velopes by comparing to IBaraffe et~a l. (2008) models of 
hot Neptunes. For planets of 10 and 20 Mq with 10% 
by mass layer of H2 and He, we fo und very good agree- 
ment between the model radii. The IBaraffe et all (2008) 
radii fall within the range of planetary radii derived from 
our model when uncertainties on the atmospheric ther- 
mal profile and energy budget in our model are taken 
into account. In other words, it is possible to choose 
values of T c g, 7, and To within the ranges described in 
§ I2.4lsuch that our model radii agree precisely with those 
from lBaraffe et al.l (|2008ft . Further, over the full range of 
atmospheric parameters considered our m odel radii de- 
viate by no more than 27% from those of Baraffc et al. 
(l2008h . 

3. TERNARY AND QUATERNARY DIAGRAMS 

In this work, we use ternary and quaternary diagrams 
to plot the relative contributions of the core, mantle, ice 
layer, and g as layer to the structu re of a differentiated 
exopla net. IValencia et al.l (|2007al ) and iZeng fc Seagerl 
(2008) also employed ternary diagrams to present the 
interior composition of terrestrial exoplanets, and pro- 
vide detailed discussions of thes e three-axis eq u ilatera l 
tria ngle diagrams. While both IValencia et al.l (|2007al) 
and IZeng fc Seagerl (|2008| ) considered three-component 
planets comprised of a core, a mantle, and water ices, 
our fiducial model also allows for a gas layer. Three- 
dimensional tetrahedron quaternary diagrams provide a 
natural extension of ternary diagrams to four-component 
systems. 

Quaternary diagrams are useful for plotting four- 
component data (w, x, y, z) that are constrained to have a 
constant sum (w + x + y + z = A = constant) . The axes 
of a quaternary diagram form a tetrahedron of height A. 
The four vertices of the diagram represent w = A, x = A, 
y = A, and z = A, while the opposing faces are surfaces 
on which w = 0, x — 0, y — 0, and z — 0, respectively. 
At each point inside the tetrahedron, the value of w is 
given by perpendicular distance to the w = face, and 
the values of the other components are defined analo- 
gously. Equilateral tetrahedrons have the property that 
the sum of the distances from any interior point to each 
of the four faces equals the height of the tetrahedron A. 
We are thus assured that w + x + y + z = A is satisfied 
at every point within the quaternary diagram. 

We use quaternary diagrams to plot all the possible 
ways a planet of mass M p and radius R p can be par- 
titioned into the four layers of our fiducial model de- 
scribed in § [2] In this case, the four-component data 
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that we are plotting in the diagram are the fractions of 
the mass of the planet in each of the four interior lay- 
ers {x core ,x mant ie-,xn 2 o,x B _/Yic) , which are constrained 
to sum to unity. The summits of the tetrahedron repre- 
sent extreme cases in which the planet is 100% iron, 100% 
silicates, 100% water ices, or 100% H/He. The face oppo- 
site the H/He summit turns out to be a ternary diagram 
for the gas- less interior compositions of the planet. 

4. RESULTS 

Our eventual aim is to draw robust conclusions about 
the composition of a low-mass exoplanet by fully explor- 
ing and quantifying the associated uncertainties. There 
is an inherent degeneracy in the planetary compositions 
that can be inferred from planet mass and radius mea- 
surements alone; for a specified planet mass, many dif- 
ferent distributions of matter within the planet interior 
layers can produce identical radii. In planet interior mod- 
els incorporating N distinct chemical layers, specifying a 
planet mass and radius each impose a constraint on the 
layer masses, leaving {N — 2) degrees of freedom in the 
allowed compositions {xi}. Further compositional uncer- 
tainties may be introduced if the planetary energy budget 
or chemical makeup are not well known and if significant 
measurement uncertainties are present in observationally 
derived parameters. 

In this work, we examine the constraints that can be 
placed on a transiting exoplanet 's interior using only 
structural models for the planet. By not employing 
planet formation arguments to impose further constrain 
the planetary compositions, our results remain largely in- 
dependent of planet formation theories. In this section, 
we apply our interior structure model to examine the pos- 
sible compositions of several example planets: CoRoT- 
7b, GJ 581d, GJ 436b, and HAT-P-llb. 

4.1. GoRoT-% 

The recent discovery of the first transiting super-Earth, 
CoRoT-7b, has us hered in a new era of exoplanet science 
ijLeger et al.ll200l . CoRoT-7b is on a 0.85359 ± 0.00005 
day period around a bright V — 11.7 G9V star. The host 
star is very active which complicates measurement of the 
transiting planet's mass and radius. By forcing the stel- 
lar radius to be i?» = (0.87±0.04) i? Q , a planetary radius 
of Rn = (1.68±0.0 9) R® is derived from the transit depth 
(|Leger et al.l 120091 ) . The planetary nature of CoRoT-7b 
has recently been confirmed by Doppler measurements 
revealing a planeta ry mass of M p — (4.8 ± 0.8) M® 
(jQueloz et al.l [20091. For the very first time, both the 
mass and radius of a super-Earth sized exoplanet have 
been measured, thereby offering the first hints about the 
interior composition of a planet in the mass range be- 
tween Earth and Neptune. 

In this section we do not consider the possibility 
that CoRoT-7b could harbor a gas layer or a signifi- 
cant water ocean. With an orbital semimajor axis of 
a = (0.0172 ± 0.00029) AU (about four stellar radii), 
CoRoT-7b is receiving an extreme amount of stellar ir- 
radiation. CoRoT-7b is most likely tidally locked, with 
a temperature of up to 2560 ± 125 K at the sub-stellar 
point assu ming an albedo of A = and no energy redis- 
tribution CLe ger et al.l 12009). Limits on the lifetime of 
a gas layer or a water ocean under such extreme radia- 
tion are discussed in § [6] We focus here on what we can 
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Fig. 1. — CoRoT-7b core mass fraction as a function of planetary 
radius. The planetary mass is (4.8±0.8) Mfg. We neglect the possi- 
ble presence of water or a gas layer and consider a two-layer planet 
comprised of a pure iron core surrounded by a Mgo.9Feo.iSiC>3 
mantle. The red, yellow and blue shaded regions denote the core 
mass fractions obtained when varying the CoRoT-7b mass within 
its 1 cr, 2 cr, and 3 a error bars, respectively. The black vertical 
lines delimit the measured radius R = (1.68 ± 0.09) Rq (dashed) 
and its 1 a error bar (dotted). 



learn about the composition of CoRoT-7b if it is a p urely 
dry, gas-less telluric planet. iValencia et al.l J2009) offer 
another point of view, considering the possibility of an 
H/He or vapor atmosphere on CoRoT-7b. 

We first examine the interior composition of COROT- 
7b under the assumption of an iron core and a mantle 
composed of silicate perovskite (Mgo.gFeo.iSiOa, approx- 
imately similar to Earth's mantle). When considering 
only two compositional layers, the measured mass and 
radius uniquely determine the two layer masses. The core 
mass fraction as a function of planet radius for CoRoT- 
7b is displayed in Figure [TJ The solid black line denotes 
the fraction of CoRoT-7b's mass in its iron core assum- 
ing the fiducial planetary mass M p = 4.8 M®, while the 
red, yellow and blue shaded regions delimit the 1 cr, 2 cr, 
and 3 a error bars on M p respectively. The measured 
planet radius and its 1 cr error bars are denoted by the 
dashed and dotted black vertical lines, respectively. An 
Earth-like composition, having 30% of its mass in an iron 
core and the remaining 70% of its mass in a silicate man- 
tle, is consistent with the measured mass and radius for 
CoRoT-7b within 1 a. If the CoRoT-7b core is not pure 
iron but also contains a light element, the core mass frac- 
tion at a specified planetary radius will be larger. Includ- 
ing 10% sulfur by mass in the iron core EOS increases 
the CoRoT-7b radius by 0.08 i?® at a core mass fraction 
of 1 (at the top of Figure QJ, while having no effect on 
the radius at a core mass fraction of (at the bottom of 
Figure [1]) . 

Our interior structure model can strengthen the ob- 
servational constraints on CoRoT-7b's mass and radius. 
With the assumption that CoRoT-7b does not have a sig- 
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nificant water or gas layer, some of the mass-radius pairs 
within MpilcTM and RpHaR (including the fiducial a 
mass-radius pair) can be ruled out because they corre- 
spond to bulk densities lower than a pure silicate planet. 
These excluded mass-radius pairs would necessitate wa- 
ter (or some other component lighter than perovskitc) 
in their composition. The fact that some la CoRoT-7b 
mass-radius pairs are excluded can be seen from Fig- 
ure [1] where the red band denoting planetary masses 
within 1 a of the measured value never fully crosses the 
Rp + Ian dotted line even at a 100% perovskite composi- 
tion. While most of this work is devoted to constraining a 
planet's interior structure from mass and radius measure- 
ments, this is an example of how limits on a planet's in- 
terior structure could be used to improve our constraints 
on a planet's mass and radius. 

The amount of iron in a exoplanetary mantle is not 
known. Earth's mantle has about 10% iron and 90% 
Mg by number fraction (Mgo.g Feo.iSiOa), but exoplan- 
ets may have varying amounts. [Elkins-Tanton fe Seagerl 
(2008a) describe an extreme example of a coreless ter- 
restrial planet in which all of the planet's iron is mixed 
in the mantle instead of sequestered in the core. To ex- 
plore the effect of varying the mantle iron fraction we 
present a ternary diagram in Figure [5] that shows the 
tradeoff between the mass of iron in the mantle com- 
pared to the mass of iron in the core. The fractions of 
the planet's mass in the Fe core, in MgSiC>3 and in FeSiC>3 
are plotted on the three axes. MgSiOa and FeSiOs are 
mixed together in the mantle as Mgi^ x Fe x Si03, where 
X is the number fraction of FcSiC^. The red, yellow and 
blue shaded regions denote interior compositions that are 
consistent with the measured planetary mass and radius 
to within 1 a, 2 a, and 3 a of the observational uncer- 
tainties respectively. All the ternary diagram except the 
high Fe corner (xf c > 0.76 — 0.86) is shaded to within 
3 a. Because FeSiC>3 and MgSiC-3 have similar densities 
(compared to the density contrast between pure Fe and 
perovskite), we have very little ability to discriminate 
the iron content of the mantle from a mass and radius 
measurement alone. Nonetheless, x contributes to the 
uncertainty in the core mass fraction of CoRoT-7b. 

4.2. GJ 581d 

We now consider the super-Earth exoplanets that are 
large and cool enough that they might retain a small 
hydrogen- helium gas layer. As an example we use 
GJ 581d, a MpSmi = 7.09 M® super-Earth with a semi- 
major axis a = 0.22 AU that is part of a multi -planet sys- 
tem around an L = 0.013 M3 dwarf star (jUdrv et al.l 
120071 iMavor et al1l2009a|) . GJ 581 is estimated to be 
8±l Gyr okfl The radius of GJ 581d has not yet been 
measured. In this section, we adopt the minimum mass 
for GJ 581d and consider two different possible planetary 
radii: R p = 1.5 and 2.0 i?®. While these values may not 
represent the properties of the true GJ 581d planet, we 
use them to illustrate how the possible presence of a gas 
layer and observational uncertainties will affect our abil- 
ity to make inferences about the interior composition of 
transiting super-Earths. 

The two putative planetary radii considered for 
GJ 581d lead to interior compositions having very differ- 
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Fig. 2. — CoRoT-7b ternary diagram. Plausible compositions for 
CoRoT-7b are shown, provided the planet has no interior water and 
no H/He layer. The fractions of the planet's mass in the Fe core, 
in MgSiC>3, and in FeSi03 are plotted on the three axes. MgSiC>3 
and FeSiOs are mixed together in the mantle as Mg 1 _ x Fe x Si03. 



ent characteristics. Ternary diagrams assuming a radius 
of Rp = 1.5 and 2.0 i?® for GJ 581d are displayed in Fig- 
ures l^a) and (b), respectively. The leftmost black curve 
in each diagram represents the locus of possible gas-less 
compositions for the stipulated mass and radius. The 
R p — 1.5 i?® planet is very dense and iron- rich; it could 
have a Mercury-like composition with 68% of its mass in 
an iron core enveloped by a silicate mantle accounting 
for the remaining 32% of the mass. In contrast, possible 
gas-less compositions for R p = 2.0 i?® are all icy planets 
with 25%-58% H 2 by mass. 

In Figure El each colored band designates a differ- 
ent gas mass fraction. For non-zero gas mass frac- 
tions (xgas), the relative contributions of the iron core, 
Mg .gFeo.iSi03 mantle, and H 2 ices to the solid inte- 
rior (inner three layers) of GJ 581d are plotted; effec- 
tively, the fraction of the planet mass in each of the solid 
layers is re-normalized by (l — ^h/Ho)- The non-zero 
width of the gas mass fraction bands in the ternary dia- 
grams is due to the uncertainty in the atmospheric P-T 
profile. Following the scheme described in § 12.41 we con- 
sider 7 = 0.1 - 10, T = 181 - 285, T cff = 73 - 93 for 
R p = 1.5i?®, and T cff = 59 - 75 for R p = 2.0i?®. 

Allowing for the presence of a gas layer on GJ 581d 
significantly increases the range of interior compositions 
that can produce the stipulated mass and radius. The 
more gas GJ 58 Id contains, the higher the average den- 
sity of the inner three layers must be to still satisfy the 
planetary mass and radius constraints. More gas results 
in an increase in the proportion of iron, as manifested in 
the ternary diagram (Figure \S§ by the fact that the gas 
mass fraction increases to the right toward the Fe ver- 
tex. An upper limit on the mass of gas that GJ 581d can 
support is reached if the planet has no H2O or silicates 
but consists solely of H/He enveloping an iron core (a 
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Fig. 3. — Ternary diagram for the solid core of GJ 581d. The 
GJ 581d minimum mass Mp = 7.09 M(g is assumed. Each diagram 
represents a different possible planetary radius: (a) R p = 1.5R®, 
and (b) R p = 2.0-Rgj. The relative contributions of the iron core, 
Mgo.9Feo.iSi03 mantle, and H2O ices to the mass of the solid 
planet bulk are plotted. The leftmost black curve represents the lo- 
cus of gas-less compositions, and gas mass fraction increases to the 
right toward the Fe vertex. The different colored bands designate 
various gas mass fractions (och/Hc) : 10 -7 (red), 10 — 6 (yellow), 
10 -5 (green), 10 -4 (blue), 10 -3 (magenta), and 10 -2 (cyan). For 
reference, the Earth's gas mass fraction is about 10 — 6 and Venus' is 
about 10 — 4 . The width of each of the colored bands is produced by 
varying the atmospheric parameters within the ranges 7 = 0.1 — 10, 
T = 181 - 285, and T cff = 73-93 (R p = 1.5R®) or T off = 59-75 
(Rp = 2.0i? e ) ■ 



composition corresponding to the iron vertex, Figure ■ 
This H/He mass upper limit occurs at 0.12% — 0.19% 
for R p = 1.5 R e and at 1.7% - 2.2% for R p = 2.0 R @ . 
These limits consider only the constraints imposed by 
the planetary mass and radius and not the lifetime of 
the gas layer to atmospheric escape. Having a gas layer 
contribute 10™ 5 of the mass of GJ 581d (for compari- 
son the Earth's atmosphere is about 10~ 6 of an Earth 
mass) increases the minimum iron core mass fraction for 
a R p = 1.5 i?© planet from the 68% gas- less value to 
74% — 78% and decreases the minimum H2O mass frac- 
tion for a R p — 2.0 i?© planet from the 25% gas-less 
value to 13% — 17%. Although a gas layer on GJ 581d 
can make at most a small contribution to the planetary 
mass, it can nonetheless have a very important effect on 
the allowed proportions of the inner three layers and on 
our ability to infer the planet's interior composition. 

So far we have only considered the inherent uncertainty 
in the composition of GJ 581d that could be inferred 
from a planetary mass and radius. In practice, observa- 
tional uncertainties also impact our ability to constrain 
the interior composition of a transiting super-Earth. For 
illustration purposes, we consider the same two putative 
GJ 581d mass-radius pairs, and assume an optimistic but 
plausible uncertainty of 5% on both the planetary mass 
and radius. Ternary diagrams plotting gas-less compo- 
sitions consistent with the planet mass and radius to 
within 1 cr, 2 cr, and 3 a are shown in Figure EJa) for 
Rp = 1.5R® and Figure HJb) for R p = 2.0R m . If com- 
positions including gas layers were included in Figure [4j 
the shaded n a regions would all be smeared out to the 
right and extended to the Fe vertex. 

Even neglecting the effect of a possible gas layer, the 
interior composition of GJ 581d is far better constrained 
for a radius of R p — (1.5 ± 5%) R® (Figure a)) than it 
is for R p = (2.0 ± 5%) i?© (Figure gfb) ) . The superior 
compositional constraints attained at the smaller plane- 
tary radius are a consequence of two effects. First, the 
R p = 1.5i?0 planet has a lower inherent compositional 
degeneracy for its fiducial (0 a) mass and radius. The 
R p — 1.5 Rq planet is dense enough that it must con- 
tain a large amount of iron, while the R p = 2.0 Rq has a 
more intermediate density and could be assembled from a 
wider range of combinations of iron, silicates, and water. 
This can be seen from the ternary diagrams (Figure |4]) in 
which the line representing the gas-less compositions for 
(M p = 7.09 M ffi , R p — 1.5i?©) is much shorter than the 
line representing the possible gas-less compositions for 
(M p = 7.09 Mq, R p = 2.0-R©). Second, the separation 
in the I a, 2 a, and 3 a contours of the ternary diagram 
are much wider in the case of R p — 2.0i?© in Figure[4jb) 
than they are for R p = 1.5i?© in Figure |Ua). The rela- 
tive uncertainty on the average planet density p is iden- 
tical (to first order) for both GJ 581d radii considered 

(ap/P w yJ(AM/M) 2 + (3AR/R) 2 = 16% J , while the 

spacings between iso-mass and radius curves on the 
ternary diagram are roughly proportional to cx Ap/p 2 
(Zcng & ScagcrJUooH). Thus, the separation in the 1 cr, 
2 cr, and 3 cr contours of the ternary diagram increases 
with decreasing planetary density. This example illus- 
trates how our ability to constrain the interior compo- 
sition of a transiting super-Earth depends not only on 
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Fig. 4. — Ternary diagram displaying plausible gas-less composi- 
tions for GJ 581d. An observational uncertainty of 5% is included 
on both the assumed mass (M p = 7.09 iWgj) and the assumed radii. 
Each diagram represents a different possible planetary radius: (a) 
Ftp = 1.5.R®, and (b) R p = 2. OR®. The red, yellow, and blue 
shaded regions denote compositions that are consistent with M p 
and R p to within 1 cr, 2 a, and 3 cr, respectively. 

the precision of our measurements, but also on the true 
mass and radius of the planet. For a given relative un- 
certainty on the average planet density, the composition 
can be best constrained for very dense planets (near the 
Fe vertex). 

4.3. GJ 436b 

GJ 436b, a hot Neptune orbit i ng a nearby M star 
(jButler et al.l 12004 iManess et all 12007ft . was the first 
known transiting intermediate-mass pla net. Since 
GJ 436b was found to transit its star by iGillon et al.l 



(2007b), substantial efforts have been made to mea- 
sure its mass and radius using photometric data 
from the Spitzer S pace Telescope (|Deming et al.l 120071 : 
Gillon et al.1 l2007al) . from the Hubble Space Telescope 
(jBean et al.1 12008ft . and from furthe r ground-based ob- 
servations (e.g. iShporer et all [20091. Here, we adopt 
values for th e prope r ties o f GJ 436b and it s hos t 
star given by iTorresI (|2007l ) and iTorres et all (|2008D . 
who employed a weighted a verage of light-curv e pa- 
rameters from ground-base d ( Gillo n et al.l l2007bf) an d 
Spitzer studies (|Deming et alJl2007l : IGillon et al.ll2007al ) : 
L* = 0.02601^" L Q , M p = 23.17 ± 0.79M e , R p = 
4.22±™X, and a = 0.02872t ~ AU 

The measured mass and radius of GJ 436b constrain 
its bulk interior composition. Allowed compositions for 
our fiducial planetary parameters (M p — 23.17A/0, R p — 
4.22i? ffi , T = 663 K, T eff = 70 K, 7 = 1) are displayed in 
Figure OH For our fiducial set of G J 436b model param- 
eters, the allowed compositions form a two-dimensional 
surface in the quaternary diagram (FigureEJa)). This il- 
lustrates the inherent compositional degeneracy originat- 
ing from an underconstrained interior model; the mea- 
sured mass and radius place only two constraints on the 
masses in each of the four interior layers. When uncer- 
tainties in the model parameters are considered, the sur- 
face of allowed compositions gains some thickness and 
spreads into a volume, weakening the constraints that 
can be placed on GJ 436b's composition (Figure HJ). Not 
all of the quaternary diagram is filled, however, even 
when both observational and model uncertainties are 
taken into account. Some interior compositions (specifi- 
cally those outside the red surfaces in Figure |6|) can thus 
be ruled out for GJ 436b. 

GJ 436b can support a range of gas mass fractions, but 
must have some gas. For our fiducial parameter choices, 
GJ 436b could be between 3.6% and 14.5% gas by mass. 
The gas mass fraction needed to produce the observed 
transit depth depends on the composition of the planet's 
solid core: water worlds with large ice layers fall near 
the minimum gas mass fraction (3.6%), while dry planets 
with iron-rich cores require up to 14.5% gas. The tradeoff 
between H/He and water contents is illustrated in Fig- 
ure[5l[b), in which the iron core and perovskite mantle are 
combined together on one axis to form a ternary diagram 
from the data presented in Figure [5][a) . In Figure Mjo) , 
the shaded wedge of allowed compositions slopes from 
near the pure H2O vertex toward increasing H/He and 
the opposite 0% water edge. Because the allowed com- 
positions span almost the entire H2O axis (from 0% to 
96.4%), the mass fraction of water on GJ 436b is poorly 
constrained by the measured mass and radius alone. 

The range of gas mass fractions that can be supported 
by GJ 436b strongly depends on the internal heat flux as 
parameterized by T c g. At higher temperatures, the gas 
layer is less dense and both the minimum and maximum 
gas mass fractions decrease, while at lower temperatures 
the gas layer is more dense and the gas mass fraction 
extremes both increase. For instance, at T ff = 113 K 
allowed gas mass fractions range from 2.3% to 11.7%, 
while at T cS = 58 K GJ 436b must be between 4.2% 
and 15.5% gas by mass. Figure [7] plots the gas mass 
fraction of GJ 436b as a function of T e g for various in- 
terior compositions (with all parameters other than T e g 
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Fig. 5. — Fiducial GJ 436b quaternary and ternary diagrams. 
The allowed compositions of GJ 436b for our fiducial choice of 
structural and atmospheric parameters (M p = 23.17M0, R p = 
4.22i? ffl , T = 663 K, T cff = 70 K, 7 = 1) are shown. In panel (a) 
we show a three-dimensional quaternary diagram plotting the frac- 
tion of the planet's mass in the iron core, Mgo.gFeo.iSiOs mantle, 
water ices, and H/He gas layer. The surface is colored according to 
the fraction of the mass of the planet found in the gas layer. Panel 
(b) displays the same data as (a) in a two-dimensional ternary di- 
agram. In panel (b) the core and mantle are combined together on 
a single axis, with the vertical distance from the upper vertex de- 
termined by the fraction of the planet's mass in the two innermost 
planet layers. The color shading denotes the relative contribution 
of the core to the total mass in the inner two layers. The width 
of the shaded wedge of allowed compositions is due to varying 
the ratio of Fe to Mgo.gFeo.iSiOs: the blue edge of the allowed 
compositions represents planets having no Fe, while the red edge 
represents planets lacking Mgo.gFeo.iSiOs. 



H/He 




Fig. 6. — GJ 436b quaternary diagram. Both model and observa- 
tional uncertainties are taken into account to determine the plausi- 
ble interior compositions of G J 436b in this diagram. Two different 
views of the same quaternary diagram are shown. The surface of 
allowed compositions for our fiducial choice of model parameters 
(M p = 23.17M3,, R p = 4.22R ffi , T = 663 K, T cff = 70 K, 7 = 1) is 
displayed in navy blue; this surface is the same as displayed in the 
quaternary diagram in Figure [5j a). To explore how uncertainties 
in model parameters weaken the constraints that can be placed 
on GJ 436b's interior composition, we vary each model parame- 
ter in turn while keeping all others fixed at their fiducial values. 
Two surfaces of the same color delimit the volume of composition 
space that is consistent with the range of values examined for each 
parameter. We consider 7 = 0.1 — 10 (cyan), To = 937 — 595 K 
(green), and T e ff = 58 — 113 K (orange). The yellow surfaces de- 
note the effect of varying the planet mass and radius within their 
1(7 observational uncertainties while maintaining all other model 
parameters at their fiducial values. Finally, the red surfaces de- 
limit the full volume of possible compositions obtained by varying 
all parameters within the ranges described above. 
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Fig. 7. — Gas mass fraction of GJ 436b as a function of T c g. 
All parameters other than T c ff (including Mp, R p , To, and 7) are 
fixed at their fiducial values. Curves for different end member com- 
positions of the solid bulk of GJ 436b below the H/He layer are 
displayed: pure iron (solid), pure perovskite Mgo.gFeo.iSiOa (dot- 
ted), pure water (short dashed), 25% iron 75% perovskite (long 
dashed), and 25% iron 50% perovskite 25% water (dot-dashed). 
The solid vertical line denotes the fiducial value of T e g = 70 K, 
while the vertical dotted lines delimit the range of T c ff values con- 
sidered (58-1130 K). 



fixed at their fiducial values). Using the formalism de- 
scribed in £12.41 we estimate T e s = 70T^2 K for a planet 
age of 6^5 Gyr; the age of the GJ 436 solar system is 
essentially unconstrained by observ ations since G J 436 
is unevolved on the main sequence ([Torres! 120071 ). Any 
constraints placed on the interior composition of GJ 436b 
will be sensitive to assumptions made about the intrinsic 
luminosity of the planet. 

Out of all the atmospheric parameters in our model, 
uncertainties in T e g have the most important effect on 
limiting the compositional constraints that can be placed 
on GJ 436b. In Figure [51 we explore the effect each 
model parameter has on the volume of allowed compo- 
sitions while keeping all other parameters fixed at their 
fiducial values. Varying T c g from 58 to 113 K expands 
the space of allowed GJ 436b compositions far more than 
varying T = 595 - 937 K or 7 = 0.1 - 10. The relative 
importance of the T e g parameter was not unexpected. 
The intrinsic luminos i ty det ermines the asymptotic be- 
havior of the Hansen! (120081 ) temperature profile in the 

radiative regime at larger optical depths (r > (T /T e g) 4 ). 
While 7 and T affect the temperature profile in the outer 
low-density low-optical-depth region of the gas layer, the 
intrinsic luminosity T c g dominates in the higher density 
inner regions of the radiative gas layer. As a result, T e g 
affects a larger component of the gas layer mass and 
exerts a larger influence on the transition to a convec- 
tive gas laye r and the entropy of the interior adiabat. 
lAdams et al.l ( 2008 ) also used the temperature profile 
from [Hansen (2008) and similarly found that T e g had 



the largest effect on their simulated planet radii. 



Observational uncertainties dominate most of the 
model uncertainties discussed above. The la observa- 
tional uncertainties on mass and radius are second only 
to the uncertainty in the planetary internal heat flux T c g 
in their effect on our ability to constrain the interior com- 
position of GJ 436b. This is evident from Figure [6] by 
comparing the yellow surfaces delimiting the volume of 
compositions obtained by varying the GJ 436b mass and 
radius within their la error bars and the orange surfaces 
denoting the effect of uncertainties in T c g . In this case, 
he range of plausible T c r would have to be constrained 
to better than about 20% of its fiducial value before the 
observational uncertainties in the planet radius would 
dominate the thickness of the volume of allowed compo- 
sitions. More theoretical work is required to model the 
cooling and internal heat flux of hot Neptunes and super- 
Earths harboring significant gas layers. Until progress is 
made in constraining T c g, improvements in the obser- 
vational uncertainties on the GJ 436b mass and radius 
will not translate into substantial improvements in our 
ability to constrain the GJ 436b interior composition. 

4.4. HAT-P-llb 

HAT-P-llb is the fir st hot Neptune t o be discov- 
ered by transit searches (jBakos et al.ll2010l). HAT-P-llb 
existence has since been confirmed by Dittmann et al.l 
(2009). Orbiting at a = 0.053Qlg-ggg| AU from a K4 
dwarf start with T eff * = 4780 ±50 K, HAT-P-llb is sim- 
ilar to GJ 436b in mass and radius: M n = 25.8 ± 2.9M® 
and R p = 4.73±0.16i? ffi (jBakos et alJl2010f >. Its host star 
is HAT-P-11 is 6.5±|;? Gyr old (jBakos et alJIMof), a s 
determined from Yale-Yonsei isochrones (jYi et al.ll2001[ ) . 
To date, HAT-P-llb and GJ 436b are the only known 
transiting hot Neptunes. 

Plausible interior compositions of HAT-P-llb are plot- 
ted in Figure [51 Figure [5] displays the surface of al- 
lowed HAT-P-llb compositions for the fiducial param- 
eter set (M p = 25.8M ffi , R p = 4.73i?©, T = 867 K, 
T e s — 66 K, 7 = 1), and also shows the effect of 
considering a range of values for each model parame- 
ter. The range of values employed for each parameter 
(7 = 0.1 - 10, T = 778 - 1227 K, T cff = 58- 86 K) was 
determined following the procedure described in § 12.41 
The lower limit on the range of T e ff values considered 
had to be truncated at 58 K to avoid having the gas-layer 
P-T profile enter an unphysical regime at high pressure 
and low temperatures (see § 12. 4|) . As for GJ 436b, uncer- 
tainties in the intrinsic luminosity of HAT-P-llb have an 
effect comparable to the 1 a observational uncertainties, 
and significantly weaken the constraints we can place on 
the planet's interior composition. 

We now attempt to compare the allowed compositions 
of HAT P-ll (Figured) a nd GJ 436b (Figu re M ■ At 
p p = 1.33 ± 0.20 gem" 3 (jBakos et alJl20iq ). HAT P- 
11b is less dense t han GJ 436b (p p = 1.69^;^ gem" 3 
iTorres et alll2008ft . HAT P-ll could thus support a more 
massive gas layer (up to 19.0%), and has a larger min- 
imum gas mass fraction (7.1%) for our fiducial choice 
of parameters. The effect of the average planet density 
on the gas layer constraints is partially mitigated by the 
higher level of stellar insolation received by HAT P-llb. 
When both 1 a observational and model uncertainties are 
taken into account, the allowed compositions for HAT-P- 



12 



Rogers et al. 



H/He 




Fig. 8. — Quaternary diagram for HAT-P-llb. Both model and 
observational uncertainties are taken into account to determine the 
plausible interior compositions of HAT-P-llb in this diagram. Two 
different views of the same quaternary diagram are shown. The 
surface of allowed compositions for our fiducial choice of model pa- 
rameters (M p = 25.8M®, R p = 4.73iJ e . T o = 867 K, T cff = 66 K, 
7 = 1) is displayed in navy blue. To explore how uncertainties in 
model parameters weaken the constraints that can be placed on 
GJ 436b's interior composition, we vary each model parameter in 
turn while keeping all others fixed at their fiducial values. Two 
surfaces of the same color delimit the volume of composition space 
that is consistent with the range of values examined for each pa- 
rameter. We consider 7 = 0.1 - 10 (cyan), T = 778 - 1227 K 
(green), and T e ff = 58 — 86 K (orange). The yellow surfaces denote 
the effect of varying the planet mass and radius within their la 
observational uncertainties while maintaining all other model pa- 
rameters at their fiducial values. Finally, the red surfaces delimit 
the full volume of possible compositions obtained by varying all 
parameters within the ranges described above. This figure is the 
HAT-P-llb analog to Figure|6]for GJ 436b. 



lib and GJ 436b overlap; it is plausible that HAT-P-llb 
and GJ 436b could both have the same proportion of 
core, mantle, water ices, and H/He gas layer. 

Our comparison between the possible interior composi- 
tions of GJ 436b and HAT-P-llb is fraught with compli- 
cations and should be interpreted with caution. Our con- 
clusions contrasting the possible interior compositions of 
GJ 436b and HAT-P-llb are dependent on the method 
used to constrain the intrinsic luminosity of the hot Nep- 
tunes (see § 12. 4|) . Our constraints on the planets' internal 
heat flux are admittedly rough and do not take into ac- 
count the influence that two different levels of stellar ir- 
radiation could have on the luminosity evolution of these 
two planets. In addition, significant scatter in the ob- 
servationally determined planetary masses and radii fur- 
ther hampers a comparative study of the transiting hot 
Neptunes' possible interior compositi ons. The GJ 436b 
radius obtained by Bean et al. (2008) using HST obser- 
vations is larger than that found from the infrared Spitzer 
light curves with a 92% formal significance, and would 
make GJ 436b less dense than HAT-P-llb. Improve- 
ments in the observational uncertainties on the mass and 
radii and in the constraints on the intrinsic luminosities 
of these two hot Neptunes are needed before we can truly 
make a robust comparison of their possible compositions. 

5. BAYESIAN INFERENCE APPLIED TO 
EXOPLANET INTERIOR STRUCTURE MODELS 

There are many model uncertainties that go into the 
interpretation of a measured mass and radius, and a ma- 
jor question is can we improve our deductions of interior 
composition from M p and R p by taking a more care- 
ful consideration of the uncertainties. So far we have 
presented our planet interior composition constraints by 
delimiting a range of compositions on a ternary or qua- 
ternary diagram. In our presentation (Figures [2j |4j |6j 
and [5]), we know that it is more likely that the exo- 
planet's true composition falls within the n a contours 
(surfaces) on the ternary (quaternary) diagram than out- 
side the contours (surfaces). We do not know quantita- 
tively, however, how likely it is that the exoplanet's true 
composition falls within the n a bounds. In this section, 
we present an approach that yields a more detailed map 
of the relative likelihoods of the interior compositions on 
the ternary (quaternary) diagram and that takes all the 
contributing sources of uncertainty into account in a for- 
mal way. 

We turn to a more technical description of precisely 
what the contours in Figures [H HI [SJ and [5] repre- 
sent, and why there is a more thorough approach. 
The n a contours (or surfaces in the case of quater- 
nary diagrams) delimit the range of compositions that 

are consistent with the measured planetary mass M p 

and radius R p (where the hats are used to distinguish 
measured values) to within their n a error bars for 
some choice of 7, To, and T e ff within the ranges de- 
scribed in § 12.41 In other words, for every composi- 
tion within the n a shadings on the ternary or qua- 
ternary diagram, there is at least one choice of the 
model parameters within the parameter space cube 

\M P - na Mp , M p + na Mp ) X \ R P - na Rp , R p + na Rp ) x 

ipimin , 'fmax) X (Tcff m j n , T e ff max ) X (Xbmm , To max ) that 
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yields a consistent solution. It is important to realize 
that the n a contours in our ternary and quaternary di- 
agrams do not represent confidence intervals. While one 
may make statements about the likelihood that the true 
planet mass and radius fall within n a of their measured 
values, our n a contours on the interior composition do 
not have a similar interpretation. This would be possible 
if only one model parameter were uncertain (for instance, 
if R p had an observational uncertainty, while M p and all 
other model inputs were known exactly). In reality, how- 
ever, there is more than one uncertainty (e.g., mass, ra- 
dius, model inputs), and a more sophisticated technique 
is needed to draw accurate composition contours that 
can be associated with a likelihood. 

Bayesian statistics provide a more rigorous approach 
to calculate how different sources of uncertainty combine 
and translate into ambiguities on the interior composi- 
tion of a planet. There are three categories of uncertain- 
ties. The first is observational uncertainties. The second 
is model uncertainties, in terms of the usually uncon- 
strained range of input parameters (see § I2.4|) . The third 
is the inherent degeneracy in interior compositions that 
yield a given mass and radius; in other words, the map- 
ping from composition to mass and radius is not one- 
to-one. Using Bayesian statistics, we can associate ev- 
ery interior composition with a "posterior likelihood" , 
a number quantifying our degree of belief that the par- 
ticular interior composition is the true interior composi- 
tion of the planet (given our limited knowledge of the 
planet, and our assumptions). The "posterior likeli- 
hood" function defined over the domain of possible in- 
terior mass distributions can then be used to draw well- 
defined contours (surfaces) in the ternary (quaternary) 
diagram for which the likelihood of the true composition 
falling within the contour can be stated. In § HI we are 
already drawing contours (surfaces) on ternary (quater- 
nary) diagrams constraining the interior compositions of 
planets; Bayesian statistics provides an alternative way 
to accomplish this. 

The foundation of Bayesian statistics is Bayes' Theo- 
rem, stated below in terms of the problem at hand (of 
inferring an exoplanets interior composition): 

p (C\D, A) (xQ (C\A) L (D\C, A) . (21) 

In the above expression, C represents the set of all 
model parameters (including interior layer mass frac- 
tions, planet mass, planet radius, 7, etc.), D repre- 
sents all the measured data we have (measure plane- 
tary mass, planetary radius, stellar mass, stellar age, 
semimajor axis etc.), and A denotes all of our assump- 
tions (spherical symmetry, differentiated planet, negligi- 
ble thermal corrections in the interior three layers, etc.). 
The function (C\A) is the prior probability of compo- 
sition/parameters C in the absence of measured data, 
given the assumptions. The priors incorporate assump- 
tions about the range of model parameters to consider. 
They may also include detailed physics; for instance, one 
could assume a planet formation theory and use it to dic- 
tate a priori which interior compositions are more likely 
than others. Next, C(D\C, A) denotes the likelihood of 
the measured data D for a given set of model parame- 
ters. Measurement uncertainties and correlations can be 
used to define the likelihood. Finally, p(C\D, A) is the 



posterior likelihood of composition/model parameters C 
given the measured data D and the assumptions A. This 
is what we hope to calculate. The proportionality con- 
stant in Equation (|21j) is set so as to ensure that the pos- 
terior likelihood p(C\D, A) is properly normalized. To 
make ternary diagram contour plots, one must marginal- 
ize (integrate) the posterior likelihood p (C\D, A) over all 
model parameters (in the set C) other than the compo- 
sitional layer mass fractions. The resulting marginalized 
posterior likelihood then represents the likelihood of a 
composition when the full range of nuisance parameter 
values is taken into account. 

We now provide two examples to illustrate how the 
Bayesian approach described in the previous paragraph 
can be applied when drawing inferences about an exo- 
planet's interior. To begin, we apply Bayesian techniques 
to the case of a solid gas-less planet having measured 
mass and radi us. We consider GJ 581d, adopting (as 
mentioned in § 14.21) the GJ 581d minimum mass M p = 
7.09 and two putative transit radii R p — 1.5R®, 
and Rp = 2. OR®. We further assume that the measured 
planet mass and radius each have associated 5% obser- 
vational uncertainties. In what follows, we reproduce 
the GJ 581d composition constraints displayed in the 
Figure 0] ternary diagrams, demonstrating how Bayesian 
statistics can be used to derive more informative and 
quantitative constraints on a transiting planet's interior 
composition. 

In this example, our assumptions A include the follow- 
ing. 

1. Our model described in §[2] is appropriate to char- 
acterize the interior structure of GJ 581d. 

2. GJ 581d does not have a significant gas layer. 

3. GJ 581d has a pure iron core. 

4. The Fe number fraction in the planet mantle is sim- 
ilar to that of the Earth (x ~ 0.1). 

5. The measurement uncertainties on the planetary 
mass and radius are Gaussian and uncorrelated. 

Our model parameters in this case are C = 
(Mp,x core ,x ma ntie), where M p is the planetary mass, 
and Xi is the mass fraction in the ith component. 
We do not explicitly include xh 2 o hi the parameters 
since it is determined by the constraint 1 = x core + 
x-mantie + £h 2 o- For a specified set of parameter val- 
ues, our interior structure model will calculate a plan- 
etary radius R p (M p ,x core ,x man ti e )- The data are the 
(putative) measured planetary mass and radius D = 

\M P ± &m p , Rp i cr p J ■ We use the measured planetary 

mass and radius with their observational uncertainties to 
define the likelihood in terms of a Gaussian joint distri- 
bution for the planetary mass and radius 

(Mpj ^core: ^mantle |-^; Aj 

= 1 c -(M p -M;) 2 /2<TM r 2 -(j?. B -.n;) 2 /2<7H r 2 

27TO- m p a- R p 

(22) 
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where R p = R p {M p ,x corei x rnant i e ). In this example, 
we take a flat prior 9 (M p , x core , x man ti e \A) oc 1, for 
which regions of composition space having equal area 
on the ternary diagram are equally likely. This prior 
is analogous to what we implicitly assumed when plot- 
ting the na contours in Figure 01 Given the assumed 
prior 9 (M p , x core , x mant ie\A), we multiply 9 and £ to ob- 
tain the posterior likelihood p(M p ,x corei x man tie\D,A) 
following Equation (|21l) . We then marginalize over the 
planetary mass M p , obtaining a posterior likelihood de- 
pending only on the interior composition, 

P (•Kcore: -^mantle \^ •> A*) 
/>oo 

, Xcore-) ^mantle |AA)dM p , (23) 

Jo 

for plotting on a ternary diagram (Figure [9]). In Figure [9] 
we show contours of constant posterior likelihood, and 
label each contour with the posterior likelihood that the 
true composition lies inside the contour (calculated by in- 
tegrating the posterior likelihood function over the area 
within the contour). When drawn in this way, the com- 
position contours in the ternary diagrams are Bayesian 
confidence regions with confidence values that should be 
interpreted as the "degree of our belief that the true 
composition of a planet falls within the contour given 
our assumptions and our observations of the planet. 

Applying a Bayesian analysis to the putative mass and 
radius measurements of GJ 581d, we extract more in- 
formative and quantitative composition constraints than 
those obtained from the non-Bayesian analysis in § 14.21 
The non-Bayesian na contours in Figure U effectively de- 
note the loci of interior compositions for the discrete 

mass-radius pairs (^M p ± nau v , Rp T n &R p ^ ■ By con- 
trast, the results of our Bayesian analysis (shown in Fig- 
ure 3]) take into account the full mass-radius relationship 
for each possible interior composition. While Figure H] 
does not give any indication of the relative plausibility of 
two different compositions within the same na contour, 
the Bayesian framework yields a posterior likelihood map 
p(x 

core i •^mantle 

\D,A) over the entire ternary diagram 
(shown by the color shading in Figure [SJ . Finally, on its 
own, Figure S] does not reveal an estimate of how likely it 
is that the true composition of the GJ 581d falls within 
its na bounds. The contours in Figure |9] are, however, 
associated with Bayesian confidence values. Comparing 
Figures [4] to [9] we see that in this case, given our assump- 
tions, the true composition of the GJ 581d should fall 
within the la contours in Figure 2] with a Bayesian con- 
fidence of roughly 75%. For the price of having to assume 
a prior 9 (C\A), Bayesian inference yields more detailed 
and quantitative constraints on a transiting exoplanet's 
composition than other analysis approaches. 

We now present a second example to demonstrate 
the effect of priors. We consider CoRoT-7b, and 
for illustration we make several different assumptions 
about the iron concentration in its mantle. As men- 
tioned in § 14. 1[ we assume that CoRoT-7b has a pure 
iron core and does not have a significant water or 
gas layer. With these restrictions, our model param- 
eters are C = (M p , x core , x Fe si0 3 ), where x Fc siO a and 
XMgSi0 3 denote the fraction of the planet's mass consist- 
ing of mantle FeSiC>3 and mantle MgSi03 respectively 




Fig. 9. — GJ 581d interior composition posterior likelihood dis- 
tribution. Only compositions without an H/He layer are con- 
sidered. An observational uncertainty of 5% is included on 
both the assumed mass (Mp = 7.09 Afg) and the assumed radii. 
Each diagram represents a different possible planetary radius: (a) 
-Rp = 1.5-R® and (b) R p = 2.Oi?0. The color shading in the 
ternary diagrams corresponds to the posterior likelihood distribu- 
tion p (xcore, ^mantle \D , A). Note that, for clarity, the two dia- 
grams have different color scales. The contours are lines of con- 
stant posterior likelihood labeled with a Bayesian confidence value 
indicating the "degree of belief given the prior assumptions that 
the true composition of the planet falls within the contour. The 
confidence value is the integral of the posterior likelihood function 
over the surface within the contour. Compare these diagrams to 
Figuresf4{a) and (b), which show composition constraints obtained 
from the non-Bayesian approach employed in § f4] under assump- 
tions identical to those used here. 

(zMgSiOa = 1 - x core - £FcSio 3 )- We proceed to calcu- 
late the interior composition posterior likelihood func- 
tion p (x core , a;FeSi0 3 \D, A) following an analogous proce- 
dure to that outlined in detail in the GJ 581d example 
above. Again assuming that the measurement uncertain- 
ties on the planetary mass and radius are Gaussian and 
uncorrelated, we define the likelihood in terms of a Gaus- 
sian joint distribution for the planetary mass and radius, 
as given in Equation (|22|) . Our choices for the priors 

9(M P i Xcore-) *^FcSi03 

| A) are described below. 
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For illustration, we consider three different choices for 
the prior 9 (M p , x core , XFeSi0 3 1^4 ) in Figure HU1 In all 
three cases, we take 9 to be independent of M p so that 
mass intervals of equal size dM p are equally likely (before 
taking into account radial velocity observations measur- 
ing CoRoT-7b's mass). 

For our first prior (Figure [TUJa)), we take a flat prior in 
which each division of mass between the Fe core, MgSi03, 
and FeSiC>3 is equally likely: 9\ (M p , x core , £FeSi0 3 \A) oc 
1. This prior is most similar to what we have implicitly 
assumed in plotting Figure [2] and corresponds to the case 
in which regions of equal area on the ternary diagram are, 
a priori, equally likely. 

Second, in Figure HOT b) we choose a prior where the 
mantle iron number fraction, is uniformly distributed 
between and 1 (with all possible values x core also 
equally likely). The Fe number fraction \ of a silicate 
material is defined as the ratio of Fe/ (Mg + Fe) by num- 
ber. For perovskite Mgi_ x Fe x Si03, x is related to the 
mass fractions XM g si0 3 and XFeSi0 3 through 

X ( x core, £FeSi0 3 ) 

_ ^FcSiQ 3 

£FeSi0 3 + (1 — X core ~ £FeSi0 3 ) (MFeSi0 3 /MMgSi0 3 ) ' 

(24) 

where /J.MgSi0 3 and At Fe Si0 3 are the molar weights of 
MgSi03 and FeSi03 respectively. Transforming from a 
uniform prior in x core and x to the variables x core and 
£FeSi0 3 we find 



core ? £FeSi0 3 |^4) oc 



9xFeSiO a 
(1 — X core ) (/iFeSi0 3 //^MgSiQ 3 ) 



(x F , 



cSi0 3 



(1 - X c 



X F eSi0 3 ) (MFeSiOs/MMgSiOs))^ 



(25) 

where the right-hand side of Equation (f25|) is the Ja- 
cobian determinant of the transformation. It is impor- 
tant to note that assuming a uniform prior probability on 
XFeSi0 3 is not the same thing as assuming a uniform prior 
probability on x ( or in other words 9\ 7^ #2)- A uniform 
prior probability on x effectively weights compositions 
having small mantle mass fractions more heavily those 
having large mantle mass fractions. This is because when 
£FeSi0 3 + XMgSi0 3 = 1 — x core is small, small increments 
in XFeSi0 3 can correspond to large changes in %■ 

For our final prior (Figure [TU]c), we take an extreme 
case for illustration, adopting a strong prior. We assume 
that CoRoT-7b's mantle has an iron fraction similar to 
that of the Earth. Specifically, we take x to have a Gaus- 
sian distribution with mean x = 0.1 and standard devi- 
ation a x — 0.1: 

9 3 (M p 5 %core-> ^FcSiO;} \A) 

(1 - X core ) (MFeSi0 3 /MMgSi0 3 ) e~ (x ~ x) /2<T * 



(xFcSiOs + (1 — Xcore — ^FcSiOs ) (/^FcSiOji/^MgSiOs)) 

(26) 

In Equation (|2"B"|) . x = X ( x core, ^FoSi0 3 )- Our third 
choice of prior has a strong effect on the posterior like- 
lihood distribution in Figure [TOfc . favoring compositions 
near the x = 0-1 line in the ternary diagram. 



The danger in the Bayesian approach described in 
this section is that one's prior assumptions will af- 
fect the compositional likelihoods, as illustrated by the 
CoRoT-7b example above. While modelers have not 
formally been using the Bayesian approach, they have 
been making critical assumptions that affect their in- 
terior composition interpretation of th e mass and ra- 
dius. For example. IZeng fe Seagerl (|2008f ) assumed a uni- 
form d istribution of allowed compositions; Valencia ct al. 
(|2007a|) excluded certain regions of the ternary diagram 
having low mantle mass fractions, but effectively pre- 
sented all rem a ining compositions as equally likely; and 
iFigueira et al.l (|2009f) used a planet formation and mi- 
gration model to predict which bulk compositions of 
GJ 436b may be more likely. The Bayesian approach 
above provides a framework in which the priors are ex- 
plicitly stated, whether they are flat or not. In this way, 
the effect on the results of choosing different prior as- 
sumptions can be quantified. 

Bayesian inference may or may not be the best ap- 
proach to interpret an exoplanet 's measured mass and 
radius. Our goal was to take into account, in a formal 
way, all the sources of uncertainty contributing to ambi- 
guities in a planet's interior composition. We have shown 
that the Bayesian approach is a way to meet this goal. 
Less formal approaches (such as that in §0} for constrain- 
ing a planet's interior composition can also be insightful, 
but one should be heedful of how multiple sources of un- 
certainty are combined together when interpreting their 
composition bounds. A problem with the Bayesian ap- 
proach is that, since the data available on any given tran- 
siting planet are limited, the priors assumed can have an 
important effect on the results. As long as the effect of 
the priors is explored and acknowledged, Bayesian statis- 
tics can help to maximize the compositional inferences we 
can draw from the limited data that we have on distant 
exoplanets. Regardless of any statistical approach taken, 
modelers must be explicit about their prior assumptions 
and about the precise significance of their compositional 
constraints. 

6. DISCUSSION 

6.1. External Constraints on Planetary Composition 

So far, we have only considered the constraints placed 
on an exoplanet 's interior composition by mass and ra- 
dius measurements alone. In this section, we discuss 
how planetary formation theories, compositional stabil- 
ity, and cosmic elemental abundances can be used to 
place additional constraints on a planet's interior com- 
position. 

Planet formation models predict that some inte- 
rior compositions are more likely to form than others. 
iValencia et al.l (|2007a[ ) considered the constraints im- 
posed by protoplanetary disk abundances, adopting the 
point that planets with large iron cores or large water 
ice layers but small silicate mantles are very unlikely. 
From the relative abunda nce of Si and F e in th e so- 
lar nebula (Si/Fe ~ 0.6), IValencia et all (j2007a[) pro- 
pose a minimum ratio of mantle to core mass. Further, 
IValencia et al.l ()2007al ) put forward that, since comets are 
dirty snowballs comprised of both volatiles (water) and 
dust, cometary delivery of water to a planet will simulta- 
neously deliver silicates to build up the planet's mantle at 
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Fig. 10. — CoRoT-7b interior composition posterior likelihood distribution. Three different choices of prior are shown. In panel (a) we 
take a flat prior on the core and FeSiC>3 mass fractions, 9\ (M p , x core , a:FeSi0 3 \A\ °c 1- ln panel (b) we adopt a uniform prior on the mantle 
iron number fraction, \ (Equation I I25II ). Finally, in panel (c) we take a strong prior on the mantle iron num ber fraction, in which \ is 
assumed to have a Gaussian distribution with mean \ = 0-1 an d standard deviation a x = 0.1 (Equation H26H V The color shading in the 
ternary diagrams corresponds to the posterior likelihood distribution p (x core , a;p e sio 3 \D, A) . Note that, for clarity, each diagram has a 
different color scale. The contours are lines of constant posterior likelihood labeled with a Bayesian confidence value indicating the "degree 
of belief given the prior assumptions that the true composition of the planet falls within the contour. The confidence value is the integral 
of the posterior likelihood function over the surface within the contour. 



the ra te of at least Si/H^O ~ 0.23 by mass. iGrasset et al.l 
(2009) choose a distribution of Mg/Si and Fe/Si molar 
ratios in the bulk compositions of their modeled plan- 
ets based on the measured abundances in a collection of 
planet-hosting stars. Our approach is to consider the full 
ternary diagram and to avoid imposing strong priors on 
the a priori relative likelihood of various interior composi- 
tions. In this way, we limit the effect of planet formation 
assumptions on the composition constraints we derive. 

Planet formation models also constrain the mass and 
composition of hydrogen and helium gas layers. In this 
paper, we have only considered a fixed ratio of H/Hc. 
Planet atmospheres may have a range different from so- 
lar, based on the atmosphere formation process. Planet 
atmospheres can originate from capture of nebular gases, 
degassing during accretion, and degassing during subse- 



quent tectonic activity. Outgassing would produce a hy- 
drogen atmosphere with negli gible helium, becau se he- 
lium is not trapped in rocks ([Heber et alj 120 07). The 
mass of the atmosphere create d from outgassing, how- 
ever, c ould have a wide range ( Elkins-Tanton fc Seagerl 
2008b). In contrast, the H/He composition of a gas layer 
captured from a nebula would reflect the composition of 
the nebula (modulated by ensuing atmospheric escape) 
and is presumably close to solar. It is conventional to 
accept that accretion of nebula gases is most important 
for massive protoplanetary cores; accretion of nebular 
gas is expected for rocky cores above 10 M m while of- 
ten ne glected for planets below 6 M s (e.g. iSelsis et al.1 
l2007al) . We show the full quaternary diagram because 
planets in the intermediate-mass range 6—10 (such 
as GJ 581d) may still accrete a significant mass of H- 
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rich gas (jAlibert et al.|[200l iRafikov 2006} and retain it 
under the right conditions. 

A natural question in exploring the interior compo- 
sition range of a hot super-Earth is whether or not 
a hot super-Earth can retain an interior water layer. 
We can set upper limits on the rate at which the low- 
mass exoplanet CoRoT-7b (M p = 4.8 ± 0.8 M m ) would 
lose H 2 0. CoRoT-7b is extremely close to its host star 
and suffering intense irradiation; the surface tempera- 
ture is 1800-2600 K at the sub-stellar point d epending on 
the p lanet Albedo and energy re distribution (iLeeer et all 
12009ft . Scaling the results of ISelsis et all (|2007aft to 
CoRoT-7b's semimajor axis and host star luminosity, we 
find a minimum w ater content lifetime of 0.07-1.0 Gyr. 
ISelsis et "all (|2007a) set upper bounds on the water mass 
loss of ocean planets around Sun-like stars, considering 
both energy-limited thermal escape driven by extreme 
UV and X-ray irradiation heating the planet exospheres, 
and non-thermal escape driven by erosion from the stel- 
lar wind. Although we do not have an upper bound 
on the water content lifetime, we do not consider the 
presence of a water layer (or of H/He) on CoRoT-7b 
because of its extreme proximity to the host star. If 
CoRoT-7b does in fact have a significant water content, 
it would be in the form of a super-fl uid H2O envelope 
with no liquid-gas i nterface (see, e . g.JLeger et al.ll2004l : 
ISelsis et al.ll2007al ). iValencia et all 42000) considered the 
possibility that CoRoT-7b might harbor a water vapor 
layer. 

Atmospheric escape is another process that is difficult 
to model but could potentially be helpful in interpret- 
ing the composition of a planet by ruling out regions 
of the quaternary diagram. It is difficult to predict at- 
mospheric escape as it depends on the detailed physical 
characteristics of the planet's atmosphere and its inter- 
action with the stellar insolation. Examples of properties 
on which atmospheric escape rates depend are the com- 
position of the atmosphere, the thermal structure of the 
atmosphere, the UV history of the host star, the den- 
sity of the stellar wind, the speed of the stellar wind, 
and the planet's intrinsic magnetic moment. In order 
to understand whether or not a planet has retained any 
hydrogen, one would have to model a specific exoplanet, 
taking into consideration the range of possibilities for 
the factors c ontrolling atmospheric escape . As an ap- 
proximation, iLecavelier Pes Etangsl (|2007l ) has consid- 
ered energy-limited atmospheric escape to estimate at- 
mosphere lifetimes. Following his approach, we estimate 
escape rates of 3 x 10 7 kgs -1 and 5 x 10 7 kgs -1 for 
GJ 43 6b and HAT - P-llb. respectively. As a second ex- 
ample, [sHiiliFall (d007H) have found that at GJ 581 A's 
current X-ray and EUV luminosity, GJ 581d should not 
currently be experiencing extreme atmospheric mass loss, 
although atmospheric erosion rates at earlier (and more 
active) stages in the GJ 581 system's lifetime are uncer- 
tain. It is fair to say that no published models conclu- 
sively detail the mass-loss history of a given Neptune or 
super-Earth exoplanet, and it is not clear for which cases 
this is even possible. 

6.2. Chemical Composition of Interior Layers 

The chemical make-up of a transiting planet's enve- 
lope, ices, mantle, and core is not known a priori. In this 
work, we have limited the chemical compositions that we 



consider for the interior layers of a planet to an H/He gas 
layer with solar composition, water ices, perovskite man- 
tle, and a predominantly iron core. We have explored 
the effect of varying the mantle iron content and of in- 
cluding a light element in the planet core (see § QJ . We 
selected our fiducial choice for the chemical compounds 
comprising the interior planetary layers in our model to 
capture the dominant materials making up the solar sys- 
tem planets. 

There are, however, several additional possibilities 
for the chemical make up of an exoplanet. For ex- 
ample, am monia ices change the EOS of Neptunes in- 
terior (e.g. iPodolak et al.l Tl995l ). Super-Earths that 
have outgass ed an extended hydrogen enve lope would 
lack helium (|Elkins-Tanton fc Seagerl [200 8b). A mas- 
sive CO gas layer in a hydrogen-poor planet would have 
a different EOS than an H/He gas layer, but because 
of its density, likely would not contribute to an ex- 
tended radius. Water-dominated "ocean" planets could 
have a va por atmosphere or even a superfluid surface 
layer ( e.g., iKuchnerl 120031: iLeger et al.l 120041: ISelsis et al.l 
l2007aft . Carbon planets will have different interior com- 
positions entirely, a s com pared to silic ate-based planets 
(jKuchner fc Seagerl 120051 : ISeager et alJl2007h . When in- 
terpreting the mass and radius for a given exoplanet, 
these other compositions should be included in the fu- 
ture. 

6.3. Planet Evolution 

We have found that uncertainties in an H/He-laden 
planet's intrinsic luminosity significantly weaken the con- 
straints that can be placed on the planet's interior com- 
position. In the case of GJ 436b and HAT-P-llb, un- 
certainties in T e ff even dominate the observational un- 
certainties on the planet masses and radii. Developing 
models to better predict a low-mass planet's intrinsic lu- 
minosity is thus an important endeavor to further our 
ability to study the interior compositions of super-Earths 
and hot Neptunes. 

Time-dependent simulations of planets as they age 
and cool could be employed to constrain the planets' 
intrinsic luminosities. Exoplanet evolution calculations 
have been pe rform ed in several p revio us studies (e.g. 
Baraffe et all 120031: iBmrows et al.l 1 20031: iChabrier et all 
20041: iFortney fc Hubbard! I2004t IBaraffe et all 120061: 
Fortnev et al.l 120071 : iBaraffe et all 120081 ). Simulations of 
the solar system giants have illustrated how complicated 
the process of predicting a planet's intrinsic luminosity 
can be. While simple models of Jupiter's evolution and 
interior structure ar e in good agre ement with the ob- 
served cooling rate ()Hubbardl 119771 ). homogeneous con- 
traction models predict intrins ic luminosities that are too 
low for Saturn ([Giiillot 2005|) and too high for Uranus 
and Neptune (|Stevensonl 1 1 982ft . 

In this work, we have subsumed an evolution calcula- 
tion by using a simple scaling relation to derive a plau- 
sible intrinsic luminosity range from a planet's mass, ra- 
dius, and age. An evolution calculation coupled with 
our planet interior model may eventually offer a more 
self-consistent approach to constrain the intrinsic lumi- 
nosities of low-mass exoplanets. The addition of a time- 
dependent cooling calculation would essentially shuffle 
our uncertainties in T c g to uncertainties in temperature- 
dependent EOSs, the planet's chemical composition, and 
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the planet's migration, geological, tidal evolution, and 
compositional histories. Poorly constrained planet ages 
lurther limit the improvements an evolution calculation 
could provide in the intrinsic luminosity constraints. 
These limitations and the added computational power 
demanded by time-dependent models motivate our use of 
an approximate phenomenological approach to constrain 
T c g. In a future paper, we plan to perform an evolution 
calculation to verify that the range of T e s values chosen 
in this work is representative of the uncertainties in a 
planet's age and history. 

6.4. Beyond Mass and Radius - Further Observational 
Constraints on Compositions 

In this paper, we have focused on the constraints 
that can be placed on a transiting exoplanet's interior 
composition using only knowledge about its mass, ra- 
dius, and stellar insolation (all properties that can be 
measured by current spectroscopic or photometric tech- 
niques). Are there other observations that can further 
restrict the range of interior compositions of a low-mass 
exoplanet? Transmission spectra during primary tran- 
sit can be used to discriminate between a planet with 
a significant hydrogen envelope an d a hydrogen-poor 
super-Earth (iMiller-Ricci et aT] [2009'). Observations will 
be extremely challenging, even with multiple transits 
(jKaltenegger fc Traubll200l iDeming et al.ll2009| ). In the 
case of close-in transiting hot Jupiters, apsidal preces- 
sion induced by tidal bulges on the planet could pro- 
duce observable changes in the transit light-curve shape, 
revealing additional information about the interior den- 
sity distribution of the plan et (through the Love num- 
ber; |RagQzzine^3yolf 2009). This idea is geared at hot 
Jupiters and it is unclear whether the effect will be signif- 
icant for terrestrial or Neptune-size planets. The poten- 
tial for an improved understanding of planetary interiors 
should provide strong motivation for the advancement of 
these observational techniques toward greater sensitivi- 
ties. 

7. CONCLUSIONS 

We have quantified how observational uncertainties, 
model uncertainties, and inherent degeneracies all con- 
tribute to the range of plausible bulk compositions for 
transiting low-mass exoplanets. We have only considered 
the constraints imposed on the composition by the mea- 
sured planetary mass, radius, and stellar insolation, and 
did not speculate on the formation history. Uncertain- 
ties in the formation, evolution, and age of the planets 
studied were encapsulated in the range of values chosen 
for the internal heat flux, albedo, 7, and mantle iron 
content. We summarize our main conclusions below. 

1. The interior compositions of CoRoT-7b, GJ 436b, 
and HAT-P-llb (the three lowest mass transiting 
planets known to date) are constrained by our in- 
terior structure model. 

• CoRoT-7b: An Earth-like composition hav- 
ing 30% of its mass in an iron core and the 
remaining 70% of its mass in a silicate man- 
tle is consistent with the measured mass and 
radius within 1 a. Large core mass fractions 
(^Fc > 0.76 — 0.86) are ruled out at the level 



of at least 3 a, but all other combinations of 
core mass fraction and iron mantle content 
(in water-less, gas-less compositions) are al- 
lowed; the planet could have no core or could 
be composed of up to 86% pure iron by mass 
and still fall within the 3 er error bars on 
M p and R p . If CoRoT-7b does not contain 
significant amounts of water or gas, some of 
the mass-radius pairs within M p ± \o~m and 
R p ± Io~r (specifically those that correspond 
to bulk densities lower than a pure silicate 
planet) can be ruled out. 

• GJ 436b: GJ 436b must have between 2.3% 
and 15.5% H/He layer by mass to produce 
the observed transit depth. These lower and 
upper limits on the GJ 436b H/He layer de- 
pend on the intrinsic luminosity of the planet. 
The water content of GJ 436b is very poorly 
constrained by the mass and radius measure- 
ments alone: GJ 436b could be completely 
dry, or could alternatively consist of up to 
96.4% water by mass. 

• HAT-P-llb: Nominally, HAT-P-llb's mea- 
sured density is lower than GJ 436b's. HAT- 
P-llb thus requires a higher minimum mass 
of gas (at least 7.1%) and can support a more 
massive gas envelope (up to 19.0% by mass). 
Comparisons between the range of plausible 
compositions for GJ 436b and HAT-P-llb are 
made difficult by the uncertain intrinsic lumi- 
nosities of these planets and by the scatter 
in the observationally determined masses and 
radii for each planet. 

2. Uncertainties in the intrinsic luminosities of low- 
mass exoplanets significantly weaken the composi- 
tional constraints that can be derived from a pair of 
mass and radius measurements. In the case of both 
GJ 436b and HAT-P-llb, the uncertainties on T cS 
dominate the observational uncertainties. Better 
constraints on T e g (possibly obtained through plan- 
etary evolution models) are required to improve our 
limits on the interior compositions of transiting hot 
Neptunes. 

3. The degree to which we can constrain the composi- 
tion of a super-Earth depends on the planet's den- 
sity. Putative planets with extreme densities (es- 
pecially those with very high densities) allow the 
tightest composition constraints (assuming similar 
observational uncertainties on M p and R p ). Denser 
planets will have smaller radii, however, making it 
more difficult to measure their transit radii with 
high precision. 

4. Quaternary diagrams provide a convenient way to 
illustrate the range of possible interior composi- 
tions for a transiting planet that harbors a signif- 
icant gas layer. They allow one to display inte- 
rior compositions consisting of four distinct com- 
ponents (in our case an iron core, silicate mantle, 
water ice layer, and H/He envelope). 

5. When constraining the interior compositions of 
transiting exoplanets, modelers must include in 
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their analysis many sources of uncertainty (model, 
observational, and inherent degeneracy). The 
Bayesian approach presented in § [5] provides a 
framework with which one can combine all the 
sources of uncertainty contributing to ambiguities 
in a planet's interior composition in a formal way. 
Given explicitly stated assumptions and the mea- 
sured planet parameters, the procedure described 
in § [5] outlines how to calculate the relative likeli- 
hood that any interior composition on the ternary 
diagram is the true composition of the planet. 
The likelihoods obtained can be strongly dependent 
on the prior assumptions made. In the Bayesian 
framework, however, the prior assumptions are ex- 
plicitly stated and so their effect can be explored 
and quantified. 



6. Allowing for the possibility of a gas layer in fu- 
ture interpretations of the mass and radius mea- 
surements of transiting super-Earths will greatly 
increase the range of possible interior compositions 
of the planet. The presence of even a low-mass gas 
layer contributing to the transit radius can signif- 
icantly alter the inferred characteristics of the un- 
derlying solid planet. Specifically, the higher the 
gas mass fraction the denser the solid planet inte- 



rior must be to compensate. 

Planetary science has come a long way toward under- 
standing planetary interiors. With Jupiter, we know that 
its bulk composition is dominated by 50%-70% hydrogen 
by mass; that the helium abundance in its atmosphere 
is somewhat below the protosolar value; that it contains 
1 M® < Mz < 39 M® of heavy elements; and that 
between and 11 M® of the heavy elements co uld be 
concentrated in a core (|Saumon fc GuillotJ[200l . For 
exoplanets, without recourse to in situ composition mea- 
surements and gravitational moment measurements from 
spacecraft flybys we will be permanently limited in what 
we can infer about the interior composition from the mea- 
sured mass and radius. Not only are the measurement 
uncertainties relatively large (2% at best on R p compared 
to ~ 0.01% for solar system planets), but models are also 
needed to map the planetary mass and radius into inte- 
rior composition and the model uncertainties are high. 
We will have to be satisfied simply knowing that we can 
quantify the wide range of exoplanet plausible interior 
compositions. 

We thank Thomas Beatty for useful comments and dis- 
cussions. We also thank Diana Valencia, whose referee 
comments helped to strengthen this paper. This work 
is supported by the Natural Sciences and Engineering 
Research Council of Canada. 
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